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Abstract 

We examine the propagation of gravitational waves in new field theory of gravi- 
tation recently proposed by Novello-De Lorenci-Luciane(NDL). This examination is 
done on a solvable case corresponding to a spherically symmetric static configura- 
tion. We show that in NDL theory the velocity of gravitational waves is lower than 
light velocity. We point out some consequences of this result and suggest a possible 
scenario for its verification. 
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1 Introduction 



The equivalence principle (EP), which states that all kinds of matter (including massless 
particles such as photons) interact in the same way with the gravitational field, gave 
to Einstein the possibility to treat gravitational phenomena as a modification of the 
spacetime geometry. 

Recently we have exploited some consequences of a field theoretical description of 
gravity with a new ingredient: we do not require that the gravitational energy satisfies 
the hypothesis of universality of the equivalence principle. This means that gravity does 
not couple to itself in the same way as all others forms of energy. The hypothesis that 
gravity-gravity interaction is the same as the interaction that gravity has with all other 
type of matter, as made in General Relativity, is not supported by direct observation, but 
by indirect consequences of the effect of gravitational field on matter. 

The NDL theory of gravity contains many of the ingredients of GR and is indistin- 
guishable as far as standard post-Newtonian tests and quadrupole emission are concerned 
(see ref. [0). The key point to distinguish the NDL theory from GR seems to be the 
velocity of gravitational waves. In order to examine the behavior of such waves we will 
study the evolution of discontinuities of the corresponding equations of motion through a 
characteristic surface S, using the Hadamard method. 

There is a general expectation concerning the possibility that we could detect gravi- 
tational waves in the next decade. Such situation is based on the great number of new 
experimental devices that many laboratories, throughout the world, are constructing. 
Many scientists are going to become involved in this enterprise. In this vein, there is a 
natural and direct way to test NDL theory by the analysis of the gravitational waves. We 
start to undertake this task in this present paper. 

In the section ^] we review the NDL theory showing its formal structure. We show 
that the gravitational radiation propagates on a effective metric distinct from the universal 
metric seen by the rest of matter and the static spherically symmetric solution. 

In the section |3] we make a brief review of Hadamard's method to study wave propa- 
gation. 

In the section |] we compare the velocity of electromagnetic waves in the vacuum and 
Gravitational waves in the NDL theory on the static spherically symmetric metric and 
show that in this theory the velocity of gravitational radiation is less than c [We use units 
of light velocity (c = 1)]. 

2 A Short Review of the NDL Theory of Gravity 
2.1 General Features 

In a previous paper ]l| we presented a modification of the standard Feynman-Deser ap- 
proach of field theoretical derivation of Einstein's general relativity, which led to a com- 
petitive gravitational theory. The main lines of NDL approach can be summarized as 
follows: 
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• Gravity is described by a symmetric second rank tensor tp^ that satisfies a non- 
linear equation of motion; 

• Matter couples to gravity in an universal way. In this interaction, the gravitational 
field appears only in the combination 7^ + tp^, inducing us to define a quantity 
9p,v = Ifiu + (Pfiw This tensor acts as an effective metric tensor of the spacetime 
as seen by matter or energy of any form except gravitational energy; 

• The self interaction of the gravitational field break the universal modification of the 
spacetime geometry. 

We define a three- index tensor F^^, which we will call the gravitational field, in 

terms of the symmetric standard variable (which will be treated as the potential) to 
describe a spin-two field, by the expression^ 

F a p^ = -^P,j,[a;P] + F[alf)]n) W 

where F a is the trace 

F a = F^-f = p >a - ¥W7 M "- (2) 

From the above definition it follows that F^^ is anti-symmetric in the first pair of indices 
and obeys the cyclic identity, that is 

Fccfiv ~\~ F^ au (3) 

Fafiv F^ ua -\- F va ^ 0. (4) 

In order to construct a non-linear theory for the gravitational field we make the hypoth- 
esis that gravity-gravity interaction is represented by a functional of the two invariants of 
the field F Q/3fl , A and B, given by: 

A = F a , u F a ^, 

B ee F^. (5) 
and the Lagrangian density is a Born-Infeld type 

^-4. (6) 

where b is a constant and U is defined by 

U = A-B. (7) 
In this way, the gravitational action is expressed as: 

S = J d 4 aV=7£, (8) 

1 We are using the anti-symmetrization symbol [ar, y] = xy — yx and the symmetrization symbol 
(x, y) = xy + yx. Note that indices are raised and lowered by the background metric 7^. The covariant 
derivative is denoted by a semicomma ';' and it is constructed with this metric. 
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where 7 is the determinant of the flat spacetime metric 7^ written in an arbitrary coor- 
dinate system. 

From the Hamilton principle we find the following equation of motion in the absence 
of material sources: 

'LuF^A = 0. (9) 

L\j represents the derivative of the Lagrangian with respect to the invariant U. 

As we remarked before, we follow the standard procedure |2], ^] and define an effective 
Riemannian metric tensor in terms of the potential ip a p, by the expression 

Qtxv = Ifxv + VV- ( 10 ) 

This relation has a deep meaning, for all forms of non-gravitational energy the net effect 
of the gravitational field is felt precisely as if gravity were nothing but a consequence of 
changing the metrical properties of the spacetime from the flatness structure to a curved 
one. The definition of the associated metric tensor is provided precisely by the above 
expression. This means that any material body (including those massless) moves along 
geodesies as if the metric tensor of spacetime were given by the equation (|10l). This is in 
agreement with the Equivalence Principle. 

2.2 The Static Spherically Symmetric Solution 

We set for the auxiliary metric?] of the background the form 

ds 2 = dt 2 - dr 2 - r 2 (d6 2 + sin 2 6d(f 2 ). (11) 

This means that all operations of raising and lowering indices are made by this metric 
7^. We remark the fact that matter^ feels a modified geometry given by Eq. (fLOf). 

A possible static spherically symmetric solution of NDL theory has only two non zero 
gravitational components of (p^ u : 

tpoo = </° = /i(r) (12) 
¥>u = <p n = -v(r). (13) 

In this case the gravitational field F a p^ reduces to 

F 100 = -- (14) 
r 

771 -F133 11/2 riK\ 

Fl22 = ^re = 2 ur -^ r (15) 

(16) 

in which a prime ' symbolizes the derivative with respect to the radial variable r. The 
unique component of the trace that remains is Fi, 

F 1 = f i'-2-. (17) 
r 

2 We note that this metric is non-observable by any form of energy. 

3 Massive or massless particles — photons, for instance — that is, any form of non-gravitational energy 
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From these we can evaluate the invariant U : 



From the equations of motion (§) we obtain 

2M r ,r, 



r 



l-(-) 4 (19) 



H = ^VbM{F(a, V2/2) + /i } (20) 



2 

in which the constant r c appearing in Eq. ( |T9| ) is given by 



9 2M 

^c 2 s — . (21) 

In Eq. (|20|) F(a, v^2/2) is the elliptic function of the first kind and the constant /xo must 
be chosen to yield the correct asymptotic limit. The quantity a is given by 



a = arcsm 



1 - 



r \ 2 



(22) 



3 Wave Propagation 



In this section we discuss wave propagation in two steps. First we give a very intuitive 
example using the electromagnetic field. We make a brief review of Hadamard's method 
and show why in non-linear cases it is better apply this treatment. Second we use this 
method on the NDL theory of gravitation 



3.1 Hadamard Treatment 

Let us consider a charged particle at rest for t < to that starts to be accelerated in t = to. 
The electromagnetic field produced by this particle was initially independent of time and 
after that it becomes a variable field. Hence, there are two regions of space-time one 
in which the field is independent of the time and one in which the field varies. These 
regions are separated by a hypersurface E. As the accelerated particle emits radiation, 
the hypersurface £ is the corresponding wavefront. 

Mathematically such a hypersurface £ is characterized by the fact that some deriva- 
tives of the electromagnetic potential are discontinuous across S. This effect comes 
from the fact that all derivatives of the potential are zero in the region E~ but not in 
E + . These disturbances are propagated by the equation of motion. 

Although the situation described above is a good image for any field disturbance, the 
equation of motion of the gravitational field is non-linear and this fact causes a great 
difference between the analysis of the propagation of a gravitational wave and of an elec- 
tromagnetic wave. Due to this, the better and simpler method to study the propagation 
of the gravitational wave is Hadamard's method, which we explain briefly in the following. 
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Let us define a hypersurface of discontinuity E by the equation: 

S : f(x a ) = 0. (23) 

The space-time is divided in two regions: S + (where if > 0) and S~ (where if < 0). The 
discontinuity of any function f(x a ) across S is defined by: 

if(x a )} = hm (f(P+) - f(p-)) (24) 

where P + is in S + and P~ is in E~. To introduce the algebraic formalism proposed by 
Hadamard we use the definition of the total differential of a function f, to state that 

[df] = (25) 



since 

[df] = lim (df(P + ) - df(P- 
It follows 

W\ = i^dx") = (26) 
but x a is a coordinate and since the coordinate system is continuous, dx a is also, hence, 

<df 



-}dx a = (27) 



dx a - 

this differential is calculated on the surface therefore dx a is parallel to it 



dx c 



oc y? ia = K a (2J 



where if >a is normal to hypersurface E. This result is valid for tensorial functions of any 
rank. We can verify this using an suitable transformation of coordinates for a particular 
case of a vector and it is straightforward extended to higher rank tensor. 



3.2 Wave Propagation in NDL Theory 

Let us review briefly the result set up in for the velocity of the gravitational wave in 
NDL theory. We represent by the symbol [J]s the discontinuity of the function J across 
the surface S. Following Hadamard we impose the discontinuity conditions for the field: 

[F^ va h = (29) 

and 

[-F/w«;a]e = ffivak\, (30) 

where k a represents the normal vector to the surface of discontinuity S. f a ^ has the 
same symmetries as the field F a ^. Taking the discontinuity of the equation of motion 
(H) we obtain: 

/%s)A*. + 2^(?7 - O^Wm = (31) 
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in which we defined 77 = F a ^^f a ^^ and ( = F^f^. In the same way, taking the discontinuity 
of the indentity 

F a f3 X ;~/ + F/3 7 x ;a + F ia x ]f3 = -5 x a F^. p + -S^Ffl.a + -5 x p F a ^ (32) 

contracted with F a/3X k' y , we obtain: 

(77 - C)k 2 + 2F af3X f^ x k a k~< + F al3X f a k p k x + FPpkpV = 0. (33) 

Working with equations ([H]) and (|33|), we arise at the following expression for the prop- 
agation: 

A; M W + AH = (34) 
in which the quantity A Ml/ is written in terms of the gravitational field: 

A" v = 2^[F^F u {aP) - F"F U ]. (35) 

Note that the discontinuities of the gravitational fields propagate in a modified geometry, 
changing the background geometry 7^ into an effective one , 

g % = ^ + A^ (36) 

which depends on the field F a p^ and also on the dynamics. This fact shows that such 
a property stems from the structural form of the Lagrangian. Thus, in the NDL theory 
the characteristic surfaces of the gravitational waves propagate on the null cone of an 
effective geometry. We remark that this geometry is distinct from that observed by all 
other forms of energy and matter, which differs from the general relativity predictions. 
This result gives a possibility to choose between these two theories just by observations 
of the gravitational waves. In the next section we will present a summary of the main 
properties of the solutions of the gravitational field in NDL theory. Then we will evaluate 
the velocity of gravitational wave in such background and compare with the GR result. 



4 Electromagnetic and Gravitational Waves in the 
Spherically Symmetric Solution 

From what we have seen, in the NDL theory all kinds of matter and non-gravitational en- 
ergy couple to gravity differently from gravity to gravity interaction. In order to compare 
the propagation of electromagnetic and gravitational waves we will proceed as follows. 

The metric that defines the structure of the spacetime in which the discontinuities of 
the electromagnetic field propagate, is provided by Eq. (|10D. We define a 4- vector such 
that 

y„<r = (37) 

where is the inverse of the metric g^, defined by 

<?^V = tfV (38) 
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Analogous to the section (|j), this 4- vector is a gradient of the hypersurface of discon- 
tinuities. The background gravitational field we are considering here, corresponds to a 
spherically symmetric and static configuration. Thus, the relation ( j37|) reduces in this 
case to 

\2 n \2 n \2 n \2 

0. (39) 



(lof {hf {hf (h) 



1 + /i 1 + v r 2 r 2 sin 2 9 

Gravitational waves propagate as null geodesic in an effective geometry gjfg given by Eq. 
d36|). The quantity is defined by Eq. (|35|) and have only two non null components: 



A 00 = -—- 40 

6 2 r 4 v ; 

A 11 = -A 00 . (41) 



Correspondingly, the effective metric is provided by 

4 /If 2 

4 /If 2 

«2 = + («) 

flS = -^Ti; (45) 
r z sm y 

Inserting these results in Eq. ( p4|) the equation of propagation of gravitational waves 
became: 

f-i 4M2 V^ 2 d 4M 'W^ ^ ( fc 3) 2 n (A a 
We can summarize this situation as: 

• Electromagnetic waves propagate in null cone of the geometry g^; 

• Gravitational waves propagate in null cone of the geometry g B ^ v . 

A question then arises: Which wave propagates faster? In order to investigate this 
problem we can run by two equivalent ways: one can either evaluate the norm of the 
vector 1^ in the geometry g c ^ u : 

(47) 

or else evaluate the norm of the vector in the geometry g^: 

KKvT. (48) 



In this vein we evaluate the quantity ([47]) and investigate the character of in the 
geometry g^ . For the solution we are considering here, results: 

11/11 -d 4M2 W (l AM2 \(l? {k)2 (/s)2 (49) 



Substituting expression (^) in the last two terms in the right hand side of the above 
equation, we obtain 

4M 2 1 \ „ ,2 I 4M 2 1 . 



The solutions of the functions /i and v are given by (|20|) and (|19|). Since we are interested 
here just in the sign of the norm of the 4-vector in the geometry determined by g^, it 
is enough to consider only the main terms. Thus, expanding the field solutions as: 

1 l + ^ + 0(r- 2 ) (51) 



1 + /i r 

' 21f +0(r- 2 ), (52) 



1 + v r 
and substituting in Eq. (|50|) result in 

2M 2 2M 2 
ll / llfo ff = — — («o) - — (ii) • (53) 

Hence, ||/|| 9cff < 0. We thus conclude that is a space-like vector in the geometry g e ^ v , 
i.e., in the NDL theory gravitational waves travel with velocity lower than light. 



5 Conclusion 

The NDL theory prediction that gravitational waves travel slower than light has some 
interesting consequences. The first one is the possibility to exist gravitational Cherenkov 
Radiation, that is the emission of gravitational radiation when a massive particle exceeds 
locally the "graviton" speed ||. Evidently, this phenomenon will put limits on how far 
from their sources cosmic rays can be found with ultra high energies, such as the ones 
searched by the AUGER Project (E > 10 20 eV) ||. How stringent these limits are will be 
addressed in future work. 

One question naturally stands up: can we obtain observational evidence that gravi- 
tational waves travel at speeds below the light velocity in the presence of gravitational 
potentials? Very high gravitational potentials can be found in the vicinities of neutron 
stars, black holes and supernova cores. However, unless we are talking about very massive 
black holes, all these gravitational sources are not sufficiently extended in size to allow 
the integration of a large enough difference in the arriving times at Earth for photons and 
gravitons supposely generated at the same place and instant. On the other hand, black 
holes of 10 6 to 10 9 solar masses, as it is believed exist at the center of galaxies such as 
M87, M51 and others [|7|], present much better conditions to delay the gravitational waves 
relatively to the electromagnetic waves. After all, the Schwarzschild radius of a 10 9 solar 
mass black hole has the size of two Astronomical Units, a distance that light takes about 
1000 seconds to cross if traveling outside the horizon. This may give us enough time for 
the potential to act on the gravitational waves, slowing them down for sufficient time, in 
order to accumulate a time delay possible to be measured with the technology of the next 
generations of gravitational wave observatories and electromagnetic telescopes in the few 
decades to come. 
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